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We report the observation of coherent heteronuclear spin dynamics driven by inter-species spin- 
spin interaction in an ultracold spinor mixture, which manifests as periodical and well correlated 
spin oscillations between two atomic species. In particular, we investigate the magnetic field de¬ 
pendence of the oscillations and find a resonance behavior which depends on both the linear and 
quadratic Zeeman effects and the spin-dependent interaction. We also demonstrate a unique knob 
for controlling the spin dynamics in the spinor mixture with species-dependent vector light shifts. 
Our finds are in agreement with theoretical simulations without any fitting parameters. 


Understanding collective spin dynamics is a problem of 
fundamental importance in modern many-body physics. 
For example, it underlies the pursuit of spintronics in 
which spin, rather than charge, is the primary carrier for 
information processing [1]. Central to the understand¬ 
ing of spin dynamics is the role of spin-spin interactions, 
which exists between identical as well as distinguishable 
spins, and their interplay with the (linear and quadratic) 
Zeeman effects. In this regard, the ultracold spinor quan¬ 
tum gas [2, 3] provides a powerful platform for investigat¬ 
ing spin dynamics due to the high controllability. So far, 
a rich variety of phenomena have been explored experi¬ 
mentally, including spin oscillations [4-11] in spinor Bose- 
Einstein condensate (BEC) and in thermal Bose gas [12], 
as well as various types of spin textures [13-16]. Very 
recently, coherent spin dynamics and giant spin oscilla¬ 
tion of a Fermi sea [17] and its relaxation [18] have been 
investigated. 

Until now, however, spin dynamics in ultracold atoms 
has been explored only in a single atomic species. Here, 
we realize a system consisting of distinguishable spin-1 
atoms of s 'Rb and 23 Na, and demonstrate well-controlled 
and long-lived coherent spin oscillations between them. 
Our study brings out several unique features of the spinor 
mixtures. (1) For collisions between two distinguishable 
spin-1 atoms, there is no exchange symmetry requirement 
and the interaction takes place over all possible total spin 
F channels[19-25]. On the other hand, only even F are 
allowed for homonuclear collisions. (2) In a single species 
spinor gas, spin dynamics is governed by the competition 
between the spin-dependent interaction energy and the 
quadratic Zeeman shift [26], while the linear Zeeman shift 
can be gauged away as magnetization is conserved. In 
the case of a spinor mixture, both linear and quadratic 
Zeeman shifts are important. (3) Because the two species 
have different electronic structures, a differential effective 
magnetic field can be generated with a spin-dependent 
ODT. Together with the external magnetic held, this can 
be used to further control the spin dynamics. 


The interaction between two heteronuclear spin-1 
bosons at a distance r can be written as [19, 20], 

U12 (r) = (a + /3fi • f 2 + 7P0) <5 (r), (1) 

where a = (g\ + 02 )/% represents the spin-independent 
interaction. The strength of spin-dependent term is given 
by/3 = (g 2 -gi)/2. The third term, 7 = (2g 0 -3g 1 +g 2 )/2, 
operates only in the total F = 0 channel via the projec¬ 
tion operator Pq. fi and f 2 label the hyperhne spin of two 
different atoms. The coupling constants gp = 2irh 2 ap/g 
are determined by the s-wave scattering lengths ap of the 
corresponding F channels and the reduced mass g, where 
h is the Planck’s constant. From previous studies [27, 28], 
we have determined (a, /?, 7) = 2'Kfx 2 ap/ gx (78.9, -2.5, 
0.06), where ae is the Bohr radius. Similar to homonu¬ 
clear spinor gases of 23 Na and 87 Rb, both /3 and 7 are 
much smaller than a [29]. The negative /3 indicates a 
ferromagnetic heteronuclear spin-spin interaction, which 
tends to align the spins of the two species along the same 
direction. 

Let us consider the collision between a 8 ' Rb atom in 
spin state |mi) and a 23 Na atom in spin state |?tt. 2 ), which 
we denote as |toi, to 2 ) in the following. Here m = ±1, 0, 
corresponding to the three Zeeman sub-levels of the / = 
1 hyperhne state. The magnetic energy associated with 
|mi,m 2 ) will be denoted as E mi,m 2 (B). The aforemen¬ 
tioned (5 and 7 terms can support several possible het¬ 
eronuclear spin changing processes as long as the popu¬ 
lation of each species and the total magnetization along 
the magnetic held B, are conserved. This is in stark 
contrast to the homonuclear spin-1 case, where only one 
spin changing process 2 |0) O |1) + | —1) is allowed. In 
this work, we focus on the following heteronuclear spin 
changing process, 

|0,-l)o |-1,0), (2) 

which is driven solely by the /3 term [30]. 

Intuitively, coherent spin dynamics of Eq. (2) can 
be understood from the interplay between the spin- 
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FIG. 1. (color online) Magnetic energy diagram for two 
heteronuclear spin changing processes driven by /? term. The 
three spin states are represented by three small arrows and 
atomic species are color coded. A E sets the energy difference 
between the two states involved. For the process in Eq. (2), 
JO, —1) | —1, 0) (blue curve), A E crosses zero at Bo = 1.69 

G, while at the same field, for the process |0,1) |1, 0) (red 

curve), A E is around 1000 Hz. In the vicinity of Bo, the small 
spin-dependent interaction energy couples the two relevant 
states (10, —1) , | — 1, 0)) and visible oscillations occur when its 
strength is comparable to A E. While the other homonuclear 
and heteronuclear processes, including 10,1) <->■ 11,0), are far 
detuned and oscillations are greatly suppressed. 
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FIG. 2. (color online) Coherent heteronuclear spin dynamics 
at B = 1.9 G. a, b, evolution of the fractional spin popula¬ 
tions of spin-1 Rb (circle) and Na (square). Red, blue and 
gray colors label the three Zeeman states | — 1), 10), and |+1), 
respectively, c, d, magnetization oscillations of Rb and Na 
show 7r-phase shift. The differences in the oscillation ampli¬ 
tudes of Rb and Na are due to their number imbalance and 
the conservation of the total magnetization, which, as shown 
in the inset, remains approximately constant in the course of 
spin dynamics. All solid lines are sinusoidal fitting to experi¬ 
mental data (see Methods). 


dependent interaction energy (/? term) and the differ¬ 
ence of total Zeeman energy between these two states: 
A E(B) = as depicted in Fig. 1. 

In analogy to a driven two-level system, when the two en¬ 
ergies are very different, the system undergoes detuned 
oscillations with large frequency but small amplitude. On 
the other hand, when the two energies are comparable, 
the system oscillates with low frequency but large am¬ 
plitude [26]. Due to the small magnitude of /3, typical 
spin-dependent mean-field energy is of the order of 10 
Hz, and as a result, visible heteronuclear spin oscilla¬ 
tions can only occur near A E = 0. As shown in Fig. 1, 
A E depends on the magnetic held B in a non-monotonic 
manner, and in particular vanishes at a Held Bq = 1.69 
G, where one expects resonant spin dynamics. This coin¬ 
cidence is a result of the slightly different Lande g-factors 
for 23 Na and s 'Rb, including contributions from both the 
linear and the quadratic Zeeman energies [29]. Near B 0 , 
homonuclear spin dynamics is greatly suppressed due to 
large quadratic Zeeman shifts and the heteronuclear spin 
dynamics for otherwise allowed spin changing processes 
are also suppressed due to large detuning. For exam¬ 
ple, the magnetic energy difference for the spin changing 
process |0,1) <-»• 11,0) is larger than 1000 Hz and will 
be substantially suppressed. Thus, working near Bq, we 
can single out the process in Eq. (2) and obtain clear 
signatures of heteronuclear spin dynamics. 

The considerations above offer only a qualitative pic¬ 
ture of the inter-species spin dynamics. Experimentally, 
we use a bulk sample consisting of an essentially pure 


23 Na BEC and a thermal gas of 8| Rb to increase the over¬ 
lap of the two clouds. This many-body system is distinc¬ 
tively different from the conventional two-level system 
since spin- and density-dependent mean-held interactions 
enter nonlinearly into the equations of motion and fur¬ 
thermore, vary in the course of spin dynamics [29]. One 
of the important consequences is the appearance of two 
magnetically tuned resonances as we shall discuss mo¬ 
mentarily. 

We produce the ultracold mixture of 23 Na and 87 Rb 
atoms in a crossed ODT with both atoms initially pre¬ 
pared in the spin state | —1) [27, 31]. To initiate the spin 
oscillations, we apply a radio-frequency (rf) Rabi pulse 
to simultaneously prepare a coherent superposition state 
with most population in | —1) and |0) for both Rb and 
Na, while populations in the 1+1) states are typically less 
than 10%. To monitor the spin dynamics, we detect the 
fractional spin population p x m = N^/N 1 for each species 
from the absorption images after various holding time. 
Here N. l m is the atom number of species i in spin state 
|m). N l = IVlj +Aq + is the total number of atoms 
of species i, with i = Na, Rb. 

Fig. 2a and 2b show typical time evolution of 
and at B = 1.9 G, respectively. The population 
in states | —1) and |0) oscillate periodically, while those 
in state |1) stays nearly constant. It is important to 
note the following features: (1) States | —1) and |0) of 
each individual species oscillate with 7r-phase shift due 
to number conservation; (2) The synchronized oscilla¬ 
tion between the two species reflects the coherent spin 
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FIG. 3. (color online) Dependences of heteronuclear spin dynamics on external magnetic field B. a, spin oscillations for Rb 
atom at B = 1.5 G (top), 1.7 G (middle) and 1.9 G (bottom). The fractional spin populations in | —1) (red) and |0) (blue) are 
monitored over long time. Away from resonance at B = 1.5 G and 1.9 G, many oscillations are observed with short period and 
small amplitude. Close to resonance at B = 1.7 G, however, oscillations are slow but with larger amplitude. The solid curves 
are for eye guiding and error bars are from statistics of several shots. The population in |+1) state is in general less than 10% 
and furthermore remains constant in the course of spin dynamics, so it is not shown here, b, magnetic energy A E as a function 
of B with zero crossing at Bo = 1.69 G. c, d, spin oscillation amplitudes c and periods d extracted from the experimental data 
and two resonances can be clearly identified in the amplitude. Solid blue lines are calculations based on many-body kinetic 
equations using experimental atomic conditions without fitting parameters. Error bars for both the amplitude and the period 
are from fitting of the oscillations and represent one standard deviation(see Methods). Mechanisms for the observed damping 
in the oscillations will be investigated in future works. 


dynamics driven by heteronuclear spin exchange interac¬ 
tion. This is even more clearly exhibited in the magneti¬ 
zation dynamics of the two species. The fractional mag¬ 
netization for each species is M 1 = (iV^ — N' L _ 1 )/N' 1 . 
The total magnetization of the system is defined as 
M = (TVff - N5_- N^)/N, where N = £. AP 
is the total number of atoms. As shown in Fig. 2c, 2d 
and the inset, A4 Na and A4 Rb are not conserved, but 
total A4 is conserved within a few percent. The coher¬ 
ent oscillations of A4 Na and AA Rb with a 7r-plrase differ¬ 
ence is a clear signature of the coherent heteronuclear 
spin exchanging process. The clean oscillations between 
the | —1) and |0) states also indicates that homonuclear 
and the other heteronuclear spin changing processes are 
greatly suppressed. 

Similar measurements of spin dynamics are performed 
for a range of magnetic fields; three examples for Rb 
are plotted in Fig. 3a. Away from Bo, fast oscillation 
with small amplitude can be observed, while very close 
to B 0 , e.g. at B = 1.7 G (middle), the oscillation is 
slow but with large amplitude. One further feature is 
worth noticing. Compare oscillations at B = 1.5 G (top) 
and 1.9 G (bottom), we note that the initial slopes of 
population change for the same spin states have opposite 
signs on different sides of Bo- This is a direct reflection 


of the sign change in A E, as depicted in Fig. 1. Sim¬ 
ilar behavior is observed for Na. These correlated spin 
oscillations for two species are well reproduced in our nu¬ 
merical simulations and furthermore are consistent with 
the ferromagnetic exchange interaction /3 < 0. 

We extract the oscillation amplitudes and periods for 
different magnetic fields and summarize the results in 
Fig. 3c and 3d. Near Bo, the system is in the inter¬ 
action dominated regime where an asymmetric double 
peak appears in the oscillation amplitude with a non¬ 
zero dip in between. This can be understood by noting 
that resonances appear when the absolute values of A E 
and spin-dependent interaction are comparable, which 
can occur on either side of B 0 , analogous to the single 
species case where A E is tuned by quadratic Zeeman 
shift [32]. However, the exact resonance positions depend 
also on homonuclear spin-dependent interactions and ini¬ 
tial conditions. The double peak behavior is, however, 
not readily distinguishable in the period where only one 
peak is observed [29]. 

To understand the observed spin dynamics quantita¬ 
tively, we model the Na condensate with the time depen¬ 
dent Gross-Pitaevskii equation [21, 30] and the thermal 
Rb cloud with the kinetic equation for the Wigner distri¬ 
bution function [12, 33, 34], The dynamics of the two 
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FIG. 4. (color online) Optical control of coherent heteronu- 
clear spin dynamics with vector light shift, a, modified de¬ 
pendence of A E for Rb on B when light induced effective 
magnetic field B ac is taken into account. Zero crossing of A E 
moves towards high magnetic fields B as B a c increases; while 
for sufficiently negative B ac , zero crossing vanishes, b, reso¬ 
nance positions as observed in the period vary with changing 
B ac , and follow the locus of the minimum of |A.E|. Here the 
external magnetic field B is applied horizontally with a 31° 
angle to the trap beams. Solid curves are for eye guiding and 
error bars are from fitting of the oscillations. B ac is calculated 
based on the measured light intensity / and p. 


species are coupled through the interaction in Eq.(l). 
Within the random phase and single mode approxima¬ 
tions [29], our simulation agrees well with the measure¬ 
ments, as shown in Fig. 3c and 3d. The simulated os¬ 
cillation frequency shows only a small kink near B 0 , in 
accordance with the experiment. 

A unique feature of the heteronuclear spin dynamics 
is its dependence on the vector light shift, which is spin- 
and species-dependent [29]. In the following, we tune 
the ellipticity of the ODT beams to further control the 
coherent spin dynamics. In the case of large detuning A, 
exceeding the excited state fine structure splitting Aps, 
the spin-dependent vector light shift is [35] 


Um{r) oc 


pm A FS 
^ A^ 




(3) 


where ojq is the energy splitting between the ground state 
and the center of the D-lines, and I(r) is the light inten¬ 


sity. The factor p characterizes the amount of circular 
polarization with p = 0 for linear and p = ±1 for pure a ± 
circular polarizations, respectively. U m can be treated as 
a “fictitious magnetic field” in the light propagation di¬ 
rection [36]. Its projection, B ac , along the quantization 
axis alters the effective magnetic field seen by the atoms. 
Due to the larger A, uq and smaller Apg for 23 Na, B ac 
for 23 Na is more than hundred times smaller than that 
for 87 Rb. For our final ODT, B ac ss 1.6 rnG for Rb and 
14 /rG for Na if p = 1. So effectively speaking, by tun¬ 
ing p, we can control the linear Zeeman energy for Rb 
and Na independently [29]. The measurements shown in 
Fig. 2 and Fig. 3 are performed with B ac essentially zero. 

Although small, B ac has a dramatic influence on the 
heteronuclear spin dynamics. For simplicity, the much 
smaller B ac for Na is ignored from now on. On the other 
hand, the much larger B ac for 8 'Rb can shift A E sig¬ 
nificantly, as illustrated in Fig. 4a. For B ac < 0, the 
zero crossing point is shifted to smaller external mag¬ 
netic fields. Eventually, for B ac < —0.2 mG the entire 
A E curve is shifted to below zero and the zero crossing 
disappears. In such cases, the spin dynamics will be es¬ 
sentially driven by Zeeman energies with a peak at the 
field of minimum |AE|. When B ac > 0, the zero cross¬ 
ing point and thus the resonance position always shifts 
to higher magnetic field. 

Experimentally, p, and hence B ac , can be tuned by 
applying the external magnetic field in the horizontal 
plane and inserting a A/4 waveplate into one of the ODT 
beams. Here p = sin(20), where 9 is the angle between 
the waveplate’s axis and the input linear polarization of 
the light. For the typical range of 9 varied in our ex¬ 
periment without causing significant heating, B ac ranges 
from —0.32 to 0.32 mG. As shown in our measurement in 
Fig. 4b, a rather small B ac can cause a significant change 
of the resonance position. For example, at B ac = 0.32 
mG, the resonance is shifted upwards by about 0.4 G. 
On the other hand, for negative B ac such that the zero 
crossing disappears, the lineshape of the oscillation be¬ 
comes much broader, as for example when B ac = —0.32 
mG, where the oscillation is always far off resonance. 

In conclusion, we have observed interaction driven co¬ 
herent spin changing dynamics between two different 
spin-1 Bose gases. Both the oscillation period and ampli¬ 
tude can be tuned over a large range with either external 
magnetic fields or, quite unique to our system, the light 
polarizations of the ODT. This latter capability is espe¬ 
cially promising because it allows sensitive and versatile 
control of the spin dynamics, as demonstrated in our ex¬ 
periment. Our system may also serve as an ideal plat¬ 
form for simulating complicated spin dynamics in solid 
state physics, such as coupled electronic and nuclear spin 
systems. 
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I. EXPERIMENTAL METHODS 

A. Preparation of the spinor mixture 

We produce the ultracold mixture of an essentially pure 23 Na BEC with l.OxlO 5 atoms and a s 'Rb thermal gas 
with about 6.3xl0 4 atoms in a crossed ODT, following the procedure described in [1], Thermal 87 Rb is used here to 
increase the overlap between two species as double BEC is immiscible [1]. The number imbalance is chosen such that 
the 8 'Rb cloud is ~100 nK above its BEC transition temperature. The final stage of the evaporation is performed in 
the presence of a 2 G magnetic field to make sure that all atoms are polarized in the | —1) spin state. The final trap 
frequencies for Rb and Na are 2-7T x (110,215,190) Hz and 27 t x (98,190,168) Hz, respectively. The average density is 
5.9xl0 13 cm -3 (6.5xl0 12 cnr -3 ) for Na (Rb). Along the vertical direction, the Thomas-Fermi radius of the Na BEC 
is 7.1/mi and the size of the thermal Rb cloud is 5.1/im. There is a differential gravitational sag of about 2.4/znr due 
to the trap frequency difference in the vertical (/-direction. 

The magnetic field is then tuned to values near Bq, which are calibrated to within ±1 mG with microwave Rabi 
spectroscopy performed on Rb. A rf pulse with Rabi frequency much larger than the quadratic Zeeman energies is 
applied afterward for spin states initialization. Three-level Rabi oscillations are observed for both species simultane¬ 
ously and we choose a pulse duration of 10/zs, which corresponds to a pulse area of 7r/3, in order to put most of the 
populations in the | —1) and |0) states. 


B. Spin population measurement 

After holding for various durations, we switch off the ODT and detect the spin population of each species with spin- 
selective Stern-Gelarch time-of-flight absorption imaging, by applying a magnetic field gradient during the expansion. 
For each experimental run, we take two images, one for each species. To obtain accurate atom numbers, our image 
system is calibrated for both Na and Rb [2, 3]. 


C. Fitting of the spin oscillations 

Spin oscillations have the same period for Na and Rb. But due to the smaller number of Rb and the the fact that 
total A4 is conserved, the fractional populations of Rb show larger oscillation amplitude than Na and offers better 
signal to noise ratio. We thus fit the Rb oscillation only with a damped sinusoidal function of amplitude A, for various 
spin components and period T 

Pi(t ) = Ci + A i e~ t l T sin(27rf/T + ip), (SI) 

where C) is an offset determined by the initial state, r is the damping constant and ip is the initial phase of the 
oscillation. Close to Bg, oscillations are typically strongly damped, in which case, only the first few periods are used 
in the fitting. 


djwang@phy.cuhk.edu.hk 
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II. ZEEMAN SHIFT AND LIGHT INDUCED MAGNETIC FIELD 

A. Zeeman shift 

The Zeeman shifts for both s 'Rb and 23 Na can be expressed as 

H z = -pF z +qF'l (S2) 

where p and q are the linear and quadratic Zeeman shifts, respectively. Both of them can be obtained from 
the power series expansion of the Breit-Rabi formula. The fine structure and nuclear Lande g-factors are gj = 
2.00233113(2.0022960) and gj = -0.0009951414 (-0.0008046108) for Rb(Na) [4], The different values of parameters 
(gj and gj) for Na and Rb renders different linear as well as quadratic Zeeman shifts (p and q ), and we have to keep 
both terms in the theoretical treatment. 

With these g-factors, the Zeeman energy difference A E ( B) (see main text) is calculated. For the spin exchanging 
process |0, —1) •o- | —1,0) which we are focusing on in this work, A E = 0 at 1.69G. 

In the presence of light induced effective magnetic field B ac , the A E = 0 point will be shifted to new values 
determined by replacing B with B + B ac for Rb. As discussed in the main text, the much smaller light induced B field 
on Na can be ignored. For Rb, pnb/B ss 702.369 (kHz • G -1 ). As B ac is typically small in our system, its contribution 
on q can also be ignored. 


B. Calculation of light induced effective B field 

As the frequency of our ODT laser is red detuned with respect to the Rb and Na D-lines with detunings A much 
larger than the fine structure splitting Aps, the optical trap potential can be approximated as [5] 

Gdipole( r ) = £4caler( r ) + U m ( I") = ^3 ^ ^ 1 + ^ fP9f m ^ ^ ^( r ) • (S3) 

Here c is the speed of light and gf is the hyperfine Lande g-factor. T and oj o are the natural linewidth and transition 
frequency of the D-lines, respectively. For the / = 1 hyperfine state, we have m = 0 and ±1. U sca \ el is from the scalar 
polarizability which causes an overall shift to all m-states. U m is from the vector polarizability and is m-dependent. 
It is only present when the circular polarization component p is not zero. Effectively, the vector potential can be 
treated as a magnetic field [6]. 

Our ODT is formed by two linear and orthogonally polarized beams propagating in the horizontal plane and crossing 
each other at an angle of 62°. To generate the effective magnetic field, we insert a A/4 waveplate into one of the ODT 
beams. The crystal axis of the waveplate is first aligned with the light polarization. Then by rotating the waveplate, 
we can control p = 8111(20), with 9 the angle of rotation. 

Experimentally, 9 has been limited to ±6° as excessive heating occurs at larger 9 , probably due to interference 
between the two ODT beams. After obtaining the light intensity from the measured laser power and beam waist, we 
calculate the effective B field as 


H ac = coscf> x U m /gfmpB, (S4) 

with (j> the angle between the laser beam and the quantization axis defined by the externally applied magnetic field. 
Here pb is the Bohr magneton. When the quantization axis is along the vertical direction, </> = 90° and thus B ac is 
essentially zero. This is the case for the measurements presented in Fig. 2 and Fig. 3 of the main text. When the 
quantization axis is in the horizontal direction, (f> = 31° and thus we can have a non-zero B ac for finite p. 


III. SOME PARAMETERS USED IN THE CALCULATION 

In this section, we give explicit expression for some parameters that enter the theoretical simulations. 


A. scattering lengths 

For the homonuclear cases, the total spin can only be F = 0 and 2. As listed in Table SI, ao and <22 are taken 
from [13] for Rb and [14] for Na. The spin-dependent interaction is determined by the scattering length difference 
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A a = a 2 — cio- We notice that uncertainties remain in Aa as reported by several groups [15-18]. The heteronucelar 
scattering lengths in Table SI are calculated based on the electronic singlet and triplet scattering lengths, and the 
molecular potentials in [11, 12]. The uncertainties of these values are on the one percent level. 



Rb-Rb 

Na-Na 

Rb-Na 

0 

II 

fc, 

= 101.8 

«N°i = 47.36 

a Rb-Na = 83.81 

F= 1 

/ 

/ 

“Rb-Na = 81-37 

F = 2 

= 100.4 

< = 52.98 

«S-Na = 76.38 


TABLE SI. Homo- and heteronuclear total spin scattering lengths used in our simulation. All the scattering lengths are given 
in a b, the Bohr radius. 

The uncertainty in these scattering lengths is one of the main causes of the small discrepancy between theory and 
experiment in Fig. 3c and Fig. 3d of the maintext. 


B. Interaction parameters 

For the homonuclear case, the two-body interaction can be expressed as 

V 12 (ri - r 2 ) = (c^ Rb + c^ Rb fi • f 2 ) <5 (r, - r 2 ). 

The interaction parameters c^ 2 Rb can be written as 


(o) 

C Na,Rb 


A fc 2 ^( 0 ) 9 _( 2 ) 

47 Th a Na,Rb + ^ q Na,Rb 

WNa.Rb 3 


( 2 ) _ 47 Th 2 QNa.Rb a Na,Rb 

^Na,Rb o 

ttlNa.Rb <3 

The parameters a, /3 and 7 can be written as 

97 r t 2 „(!) ,„( 2 ) 

znn “Rb-Na T M Rb-Na 


/? = 


g 2 

9^2 J 1 ) ,„( 2 ) 

znn URb-Na ~r “-Rb-Na 

Li 2 

27 tH 2 2ai b _ N — 30,,, 


,( 2 ) 


7 = 


Rb-Na ' Rb—Na 




where /.i = TON a ttiRb/(rtiNa + ffiRb) is the reduced mass. 


(55) 

(5 6 ) 

(57) 

(SB) 

(S9) 

(S10) 


IV. THEORETICAL CALCULATIONS 

In the main text, a simple picture based on two-level system is introduced. Such a picture is helpful for establishing 
an intuitive understanding of the spin dynamics, but cannot catch all the many-body effects. For a complete descrip¬ 
tion, we have developed a mean-field model. Here both the two-body and the many-body theories are presented. 


A. Two-body theory 

Let’s focus on the spin changing process [0, —1) O | —1,0), and ignore influences of all other processes. As M = 
mi + m 2 = — 1 is conserved, the above spin changing process involves only total spin F = 1 and 2 channels. The 
spin-dependent interaction Hamiltonian can then be expressed as 

Hi = ^2,g F | F, M) (F, M\ = g 2 |2, -1) (2, -1| + 9l |1, -1) (1, -1|. 

F 


(sii) 
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Evaluate the Clebsch-Gordan coefficient between \mi,m 2 ) and \F,M), we have 

| 2 ,- 1 > = -^(| 0 ,- 1 > + |- 1 , 0 » 

|1,-1> = -L(|0,-1>- 1-1,0)). 

Thus in the basis of |0, — 1) and | —1,0), H\ can be expressed as 

92 + 9i 92- 9i 

92- gi 92 + 91 J ’ 

where (g 2 — g\)/2 is exactly the spin-dependent /? term. The Zeeman energy can be written as 



H z 


E 1 (B) 0 

0 E 2 {B) ’ 


(512) 

(513) 


(S14) 


(S15) 


and the total Hamiltonian is 


H = H l +H z 


{92 + gi )/2 + Ei(B) (g 2 -gi)/2 

{92 — 91 )/^ (g 2 + gi)/2 + E 2 (B) 


(S16) 


The generalized Rabi frequency of the two-particle system is given by 

ft = V(E 1 (B)-E 2 (B)Y + (g 2 -g 1 )yH = y/AE(B)* + (g 2 - g^/H. (S17) 

At B 0 = 1.69 G, where A E(B) = 0, this model predicts resonant oscillation with the longest period and largest 
amplitude. This result disagrees with our observation apparently. 

On either side of Bq , we should have detuned oscillations with shorter periods and smaller amplitudes. When A B 
is much larger than g 2 — gi, the amplitude becomes too small to be observed. Our measurements in this region agree 
with this prediction qualitatively. The oscillation period can be well approximated as I/O. 


B. Many-body theory 


In this part, we outline the basic kinetic approach that we employ to interpret the experimental data. For s 'Rb 
thermal cloud, the evolution of spin polarization as well as spatial distribution can be described, within semiclassical 
approximation, with the so-called Wigner function 

g i:j (R,p-,t) = J d 3 re~ ipr (Vj(R+ |,i)^ , (S18) 

where V’*( r ) is the Heisenberg annihilation operator for Rb atom with spin a and at position r. In our calculation, 
all three spin components of the F = 1 manifold of the 8l Rb are taken into account and as a result, g l3 (R„ p: t) is a 
3x3 matrix. In the following, we shall label the Zeeman levels by i,j = — 1,0,1 for both Na and Rb. The equation 
that governs the time-dependences of the Wigner function g i3 { R, p; t) can be written in the following general form [7] 

dg(R _ P;t) + . V < 7 (R, p; t) - ^{V R I/ Rb , V p 5 (R, p; t)} - Ug{ R, p; t), U Kh ] = I Rb . (S19) 

ut nr Rb z n 

where / Rb is the collision integral that describes the effects of interactions that are not captured in the effective 
potential E/ Rb (R, p;t), which is in general a (matrix) function of R, p and t. In our calculation, t/ Rb (R, p;£) is 
obtained within the random phase approximation (RPA), which we generalize to the case of Bose-Bose mixtures. The 
anti-commutator { A , B} and the commutator [A, B] refer to quantities in spin space. 

For the Na condensate, we shall make use of the mean-field Gross-Pitaevskii approximation. Let the field operator 
for the Na condensate to be <j>i( r, t), then we can replace it by its expectation value t) —> (<t>i( r, t)), which we will 

denote simply as bi(r,t). To find the expression for G Rb (R, p;t) and the time dependences for bi(r,t), we shall start 
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with the full Hamiltonian, which can be written in three parts: 

K 2 X / 2 


H Rb = J d 3 rjv’J 


^Na = / d 3 r 


'3„ ) B 


2 TORb 

n 2 v 2 


Sij — PR b F z + qRbF 2 + y Rb 


Sij PN a. F z ' </\ a F z -\- V^sj a 

t 




fij + tU] 


So) 


s 2 ) 


^f-SaSjk + ^F u ■ Fj k 


SO) 


S 2 ) 


C -fS u 6 jk + ^F a ■ F jk 


il’kin 


( _1 \i—j 

•ffRb-Na = / d 3 r <{ aipjipifljtfrj + Pil>]Fuij)i ■ 4>\F jk <j) k + 7-— - V’I'lM-s'/’-j 


(520) 

(521) 

(522) 


Here, p Rbj Na and g Rb) Na are the linear and quadratic Zeeman energy for Rb and Na atoms. VR bj Na are the confining 
harmonic trapping potential for the Rb and Na components, which we assume to be spin-independent. are the 

standard interaction parameters for the Rb and Na spinor gases [8]. In the case of Bose-Bose mixture, inter-species 
interactions are described by three independent interaction parameters a,/3 and 7 [9]. Their explicit expressions in 
terms of scattering lengths are given in Sec.III. 

There are three different contributions to I7 Rb (R, P; t). The first comes from the single particle term 

t4b W = ^Rb (R) - PRbF z + qR h F 2 . (S23) 

The second term comes from the intra-species interactions. Here one uses the RPA approximation and obtain the 
following expression [7] 

[/^(R; t) = Cg b (Tm Rb + n Rb ) + c^ b Tr(F?r Rb ) • F + c^ b Tr(F?r Rb • F). (S24) 

where F = (F x , F y , F z ) are the spin-1 operators and n Rb is defined to be 

nRb,ij(R,t) = ^j(R^)^i(*M)) ■ (S25) 

The last contribution comes from the inter-species interactions and can be written as, within GP mean-field theory 

U™ (R; t) = aTr(riNa) + /3Tr(Fn Na ) ■ F + 7«*, (S26) 

where we have defined 


and 


nNa,jj(R,i) = bj(R,t)bi(R,t) 



bUb-i ~bUbo bUh 
-&0&-1 b* 0 b 0 -b*^ 

bib., -blb 0 bib ! 


Similarly, we can write the equation of motion for the condensate order parameter b' 


[bU^bl). 


i 4t b = 


H 2 X7 2 

2m Na 


- Pn&Fz + qu a F 2 + F Na + c^ a Tr(?r Na ) + c^Fft) • F 


aTr(n Rb )& + /3Tr(n Rb F) • F b + jU^b. 


Here U^ is defined to be 


IStp — 


1 

3 


(V’liV’-i) -(V’-iV’o) (^liV’i) ’ 
-<V4^-i> (V’oV’o) -(ipoipi) 
(ipi^-i) -{'•Pii’o) (V’lV’i) . 


(S27) 


(S28) 


(S29) 


(S30) 


Now, equations (S19) and (S29) furnish a complete description of the Bose mixture system with one component being 
thermal and the other condensate. In general, they are not very easy to solve. As a result, we shall make use of the 
single-mode approximation in the following analysis. For the Rb thermal cloud, we write 


nRb,ij (R, p; t) = Z exp 


p 2 /2m,R b + F Rb (R) 

fc R T 


aij(t). 


(S31) 
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Namely, the R and p dependences of n R \,.jj is given by thermal distribution, independent of time, while its spin 
dependence is given by a t j (t ). Here 


/> R /' exp 

p 2 /2m Rh + V Rh (R)' 

/ / (2ttA)3 6 P 

k bT 


is the classical partition function. For Na condensate, we write 

bi(R-,t) = >/Tr(n N a(R))CiW = \/ra c (R)Ci(i), 


(S32) 


(S33) 


where the total condensate density of Na, n c (R), is independent of time. The spin part of the condensate wave 
function is given by Q(t) which carries all the time-dependences. The matrix a and the spinor Q : satisfy the following 
conditions: 


Tr cr = 1, ^CCi = l- 

i 


(S34) 


We shall use a Thomas-Fermi form for the condense density distribution in the harmonic trap (see Sec.Ill), 


n c (R) 


n c { 0) 


x 2 

M 


(y + vo ) 2 _ 

Rl RjJ ’ 


(S35) 


where y o = g/oj^ h is due to gravitational field, and R x , R y and R z are the Thomas-Fermi radius of the condensate 
cloud and n c (0) is the central density. 

With these preparations, we can substitute eqn.(S31) into eqn.(S19) and integrate over R and p on both sides. Let 
us neglect the collision integral I Rb , which contributes to the damping of the oscillations and will not be discussed in 
this work. The second and third terms in the left hand of eqn.(S19) vanishes upon integration over R and p. As a 
result, we find 


where we have defined 


da 

di 


i 

h 


[cr, Mtg], 


(S36) 


Mtg = -Pm>F z + q Rb F z + c™ b nTr(Fa) ■ F + c™ b nTr(Fcr • F) + /3n tc 


AW 

NRb 


Tr(Fr) • F + ^n u 


A~Na 

NRb 


U c , 


where AAb.Na. are the numbers of atoms for Rb and Na. n and n tc are defined by 


n = 


N Rb 


1 


n tc = 


d 3 R[Tr( ? r Rb )] 2 , 

d 3 RTr (riRb) n c (R). 


VNRbNNa 

We have also defined the matrix r in terms of its matrix elements 


Tij = (*(i, Trr = 1, 

and the matrix 

C-iC -1 -C-iCo C-iCi 
-CoC-1 CoCo -Co Ci 

crc-i -Ci*Co cr Ci 



(537) 

(538) 

(539) 

(540) 

(541) 


The equation for the condensate order parameter b can be written in terms of r, in a similar form as for a. By 
substitute eqn.(S33) into eqn.(S29) and integrate over R, we find 

dt T = h^ T ' ^ BEC ]’ (S42) 

with 


Al~Rb 

A^Na 


Tr(Fcr) • F + 7 n u 


N Rb 




U a , 


M B ec = Pn&Fz + Q'Na-Pf + c^ a n c Tr(Fr) • F + /3n t , 


(S43) 
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FIG. SI. Simulated heteronuclear spin oscillations for Rb (a) and Na (b) based on our theoretical model at three different 
external magnetic fields B = 1.55 G (top), 1.7 G (middle) and 1.8 G (bottom). The 7r-phase differences between oscillations of 
| —1) (red curve) and |0) (blue curve) states of the same species and the same spin state of different species, are well reproduced 
in our simulations. 


where U a is defined to be 


and 



CT—1,—1 ~<7-1,0 a -l,l 

— <70,-1 <7o,0 —00,1 

01,-1 — 01,0 01,1 


n c = 


1 

^Na 


/ 


d 3 R?i c (R) 2 . 


(S44) 


(S45) 


Equations (S36) and (S42) describe how the spin evolves under the influences of the effective external “fields”, Mtg 
and Mbec- 


C. Numerical simulations 

In this subsection, we discuss in some more detail the numerical calculations of the spin dynamics in a mixture of 
spinor Bose gases, based on the many-body theory outlined in Sec.IVB. We shall focus on two aspects of the spin 
dynamics in this section: (1) the period, amplitude and phase of the oscillation and (2) the double peak feature of 
spin resonance. 

In Fig.SI, we present the numerical simulations for three different magnetic fields B = 1.55G, 1.7G and 1.8G. To 
the left of the resonance, B = 1.55 G, the oscillation is far off resonance, which results in large frequency oscillation 
with, however, smaller amplitude. We also note that the phase of the oscillation are consistent with the experimental 
observation. In Fig.SI, we also show that the populations in the remaining states tof = 1 stay almost a constant 
throughout the simulation. Similarly, for B = 1.8 G, to the right of the resonance, the oscillation frequency is high 
while the amplitude is small. The initial oscillation phase of the tof = 0 is opposite to that at B = 1.55 G and is 
consistent with experiment. Close to resonance at B = 1.70 G, the amplitude of the oscillation is very large while the 
frequency is very low. All the above are qualitatively consistent with two-particle picture. 

A qualitatively new feature emerges close to the resonance, where the single resonance peak predicted by the two- 
body theory split into two. This feature depends on the initial state of the spinor mixture and furthermore, depends 
on the relative population of the two species. In Fig.S2, we show the amplitude and frequency of the spin mixing for 
an equal mixture of Rb and Na atoms. Clear double peak structure are observable in both amplitude and frequency. 
While in our experiment, the double peak is only apparent in the amplitude due to the number imbalance. 

Finally, in our numerically simulation, the effects of damping have not been included. However, one feature that is 
common to our experimental observation is that the damping has the strongest effect close to resonance. This can be 
explained by noticing that the effects of spin-spin interaction is the strongest, as compared with the single particle 
Zeeman energy in this region. As a result, one expects significant damping there. Other sources of damping may 





































FIG. S2. Simulated spin oscillation amplitudes (a) and periods (b) with the magnetic field near Bo with equal numbers of Na 
and Rb atoms. Na and Rb have the same oscillation period and amplitude here. Double peaks are clearly shown in this case. 


come from the residual magnetic field gradient in our experiment, which tends to drive the sample into incoherent 
mixture through spin diffusion. This later effect is very important in our system also since the Na condensate cloud 
is much larger than the Rb cloud and spin density gradient are very large especially close to the interface of the two 
cloud. We intend to study these effects more carefully in a separate publication. 
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